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In this paper we present calculations of the SO (5) quantum rotor theory of high-T c superconduc- 
tivity in Zeeman magnetic field. We use the spherical approach for five-component quantum rotors 
in three-dimensional lattice to obtain formulas for critical lines, free energy, entropy and specific 
heat and present temperature dependences of these quantities for different values of magnetic field. 
Our results are in qualitative agreement with relevant experiments on high-T c cuprates. 



I. INTRODUCTION 

The SO (5) theory of high-T c superconductivity, unify- 
ing antifcrromagnetism and superconductivity was pro- 
posed to combine SO (3) antiferromagnetic (AF) stag- 
gered magnetization with superconducting (SC) two- 
component-. U(l) real paring field into a new order 
parameter.^ The smallest symmetry group meeting this 
requirements is the SO (5) group. It contains SO (3) group 
of spin rotation and U(l) gauge group of "charge rota- 
tion" as subgroups along with additional, the so-called 
'V operators rotating AF to SC state and vice versa. 
In the SO(5) theory ordered phases arise once SO(5) 
symmetry is spontaneously broken: with SO (3) symme- 
try breaking AF phase arises, while for U(l) symmetry 
breaking SC phase appears. Consequently, AF and SC 
order parameters are grouped in five component vector 
called "superspin" , which direction is related to the com- 
petition between AF and SC states and the kinetic energy 
of the system is that of a SO(5) quantum rigid rotor. The 
low-energy dynamics is determined in terms of the Gold- 
stone modes and their interactions specified by the SO (5) 
symmetry. 

Although, the SO (5) theory was originally proposed 
in the context of an effective quantum non-linear a 
model (NLQcrM) description of the SO (5) rotor model 
its predk±janjs_ were tested also within microscopic 
modelsMMral Calculations based on NLQcrM showed 
that the features of the phase diagram deduced from 
the SO(5) theory agree qualitatively with the global 
topology of the -pbserved phase diagrams of high-T c 
superconductors.u Magnetic correlation functions within 
the SO(5) theory were investigated showing that the 
theory predicts a scenario for the evolution of mag- 
netic behaviour, which is qualitatively consistent with 
experiments.^ Furthermore, the study of the quantum 
critical point scenario within the concept of the SO (5) 
group showed that the scaling of the contribution to the 
electrical resistivity due to spin fluctuations, displayed a 
linear resistivity dependence on temperature for increas- 
ing quantum fluctuations, which is a hallmaijk-i example 
of anomalous properties of cuprate materialsJi3 Finally, 
thermodynamic properties of the SO(5) model were stud- 
ied, where entropy and specific heat were calculated and 



compared with experimental findings. 

Many experimentally observed properties of high-T c 
cuprates show strong dependence on the magnetic field, 
e.g. entropy, specific heat, magnetic susceptibility, elec- 
trical resistivity, etc. Therefore, a proper theory of high- 
Tc superconductors must be able to explain the magnetic 
properties of these materials. Since, SO (5) model in the 
presence of finite chemical potential and a finite Zeeman 
magnetic field has the exact SU(2)xU(l) symmetry, the 
Zeeman magnetic field can be introducecL-and the SO (5) 
theory can be tested in any doping level.ta 

The aim of this paper is to study influence of the Zee- 
man magnetic field within the SO (5) theory. Obtained 
results (e.g. specific heat) can be then (qualitatively) 
compared with the relevant experiments and test the ba- 
sic principles of the SO (5) theory. 

The outline of the reminder of the paper is as follows. 
In Section II we introduce the quantum SO (5) Hamil- 
tonian in the Zeeman magnetic field and find the corre- 
sponding Lagrangian of NLQcrM. In Section III we estab- 
lish the phase diagram of the system in an applied mag- 
netic field. Section IV is devoted to the study of magnetic 
dependence of various thermodynamic functions: free en- 
ergy, entropy and the specific heat. Finally, in Section 
V we summarize the conclusions to be drawn from our 
work. 



II. HAMILTONIAN AND THE EFFECTIVE 
LAGRANGIAN 

We start from the low-energy Hamiltonian of super- 
spins rii placed in the discrete 3D s.c. lattice (3DSC) in 
the Zeeman magnetic field B along y axis. The sites are 
numbered by indices i and j running from 1 to N - the 
total number of sites. The superspin components, labeled 
by fi and v — 1,...,5 refer to AF (n^ F = [n%, 713, ruL) 
and SC (nf c = (ni,rt 5 )-) order parameters, respec- 
tively. The superspin components are mutually commut- 
ing variables!] and their values are restricted by the rigid 
rotator constraint n? = 1. The SO (5) Hamiltonian 

H = ^EL^nr-E^.^ 

i fl<u i<j 
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is 



(1) with £ being the Lagrangian: 



consists of three parts: the kinetic energy of the rotors 
(where u is an analogue of moment of inertia), inter- 
site interaction energy (with Jy being the stiffness in the 
charge and spin channel) and SO (5) symmetry breaking 
part (including the Zeeman magnetic field acting in AF 
sector). The quantities = rh ili ip V i — n V i'p i a. are genera- 
tors of the Lie SO(5) algebra (related to the total charge, 
spin and so-called l 7r' operators!!!) and p^i are linear mo- 
menta given by: 
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Furthermore, we assume that Jij = J(Ri — Rj) is nonva- 
nishing for the nearest neighbours and its Fourier trans- 
form 
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(3) 



is simply J q = Je(k), where 



e q = cos q x + cos q y + cos q z (4) 

is the structure factor for the 3D s.c. latticed! 

The last three parts of the Hamiltonian provide SO (5) 
symmetry breaking terms. In the result of their inter- 
play, the system favours either the 'easy plane' in the 
SC space (nx,n 5 ), or an 'easy sphere' in the AF space 
(712,713,714). Two of the three terms influence directly 
AF order parameter: 
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where w is the anisotropy constant, B is the Zeeman 
magnetic field and L 24 is y component of the spin vec- 
tor. Positive values of w and B favour the AF state. 
The remaining term acts on SC sector and contains total 
charge operator Lj 5 , whose expectation value yields the 
doping concentration and the chemical potential [i (mea- 
sured from half-filling) , which positive value fovours the 
SC state. 

We express the partition function Z = Trer^ H us- 
ing the functional integral in the Matsubara "imaginary 
time'V formulation^ (0 < r < l/k B T = f3, with T being 
the temperature). Explicitly: 
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The problem,can be solved exactly in terms of the spher- 
ical model.t 4 ] We note that the superspin rigidity con- 
straint (nf = 1) implies that a weaker condition also 
holds, namely: 



1 N 

i=l 
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= 1. 



(8) 



Therefore, the superspin components (r) satisfying 
the quantum periodic boundary condition (/3) = (0) 
will be treated as continuous variables, i.e., —00 < 
rii (r) < 00, but constrained on average (due to Eq. (§])). 
The constraint can be implemented using Dirac-<5 func- 
tion 5 (x) = J 00 [d\/2Tr]e lXx , which introduces the La- 
grange multiplier A (r) adding an additional quadratic 
term (in fields) to the Lagrangian (Q). Consequently, 
the partition function reads: 



Z = 
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where the function <f> (A) is defined as: 
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In the thermodynamic limit (N — > 00), the method of 
steepest descent is exact and the saddle point A (r) = Ao 
satisfies the condition: 



WW 



SX(r) 



= 0. 



(11) 



A=A 



At the antiferromagnetic and superconducting phase 
transition boundaries the corresponding order suscepti- 
bilities become infinite (see, Ref. ^), which implies for 
the Lagrange multipliers: 
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(12) 



for AF and SC critical lines, respectively. 
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III. PHASE DIAGRAM 

Providing the spherical condition (p|) with values of the 
Lagrange multipliers (|l^) one can finally arrive at the 
expression for the critical lines separating AF (or SC) 
and QD (quantum disordered) states: 
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where: 

Mt) = 
B- (?) 
B + (?) 



2 A 



c(0 = >/^ 
^ - 2 M , £_(0 = Vz^zz 

(14) 

and A = Ag F (Xq C ) for AF (SC) line. For convenience, 
in order to perform momentum integration over the 3D 
Brillouin zone, we have introduced the density of states 



p(0 = ^J2 6 ^-<^ 

q 



(15) 



(for explicit formula, see Ref. ||). 

Regions of AF and SC phases are separated by the first 
order transition line (for /i = fi c ) given by the condition 
of equality of free energies of both states, which reads: 

4 ' 



\ au — \ 
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It is instructive to present the phase diagram as a func- 
tion of physically measured quantity-the charge concen- 
tration x instead of the chemical potential. In the SO(5) 
theory, the charge concentration can be deduced from the 
free energy: 



1 #(Ap) 
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(17) 



The temperature-charge concentration phase diagram 
is depicted in Fig. 1. Here, AF to SC transition line 
splits into a region of constant chemical potential, where 
AF and SC states coexist (mixed region M). This be- 
haviour can be explained by a two-phase mixture with 
different densities at first-order phase transition. In this 
case the system globally phase-separates in two different 
spatial regions with different charge densities, but the 
same free energy. As a result, the added charges only 
change the proportion of the mixture of the two phases, 
bun not the free energy, which implies infinite compress- 
ibility (defined as dx/dfj,). 



IV. THERMODYNAMIC FUNCTIONS 
A. Free energy 

The free energy is defined as / = - (f3 N)~ 1 In Z = 
((3) 1 4> (Ao). Using the formula (|Io|), it can be explicitly 
written: 

/ = -A + Fx (A) + Fx + Fx {B+) + 

+ Fx (£>_) + Fx (D+) , (18) 

where the function: 



Fx (H) = I jT p (?) d? In {2 sinh (?) j 



(19) 



and A (?), B-(£), B+ (?), L>_ (?) and D+ (?) are defined 
by formulas (|l4]). 



B. Entropy 



The entropy is defined as S = /cs/3 2 |^. Using the 



formula (llq) we obtain: 



S (/3) = F 2 [A) + F 2 + F 2 (B+) + 

+ F 2 (D-) + F 2 (D+), (20) 

where the function: 



F 2 (E) = ^ / r t<u!' [ ./EiiH-'l. 
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C. Specific heat 



The specific heat at constant volume is defined: 
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The temperature dependence of the specific heat in var- 
ious magnetic field is depicted in Fig. 2. For low tem- 
peratures C (T) is independent on B and can be approx- 
imated by C (T) ~ T 3 . For higher temperatures the 
dependence become roughly linear with higher values of 
C (T) for higher magnetic fields. The critical tempera- 
ture is depressed for higher fields, which is in agreement 
with experimental findings. Only, the size of finite jump 
occurring at the critical point, seems to be independent 
on B (for Y-123 jump of specific heat in critical temper- 
ature is absent because of flux lattice melting, see Ref. 
[l5|,|l6]). Consequently, Zeeman magnetic field does not 
influence the value of the critical exponent a — 0. In 
high temperatures C (T) dependence saturates. 
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V. FINAL REMARKS 

In conclusion, we have considered the influence of the 
Zeeman magnetic field in the SO (5) theory of high-T c su- 
perconductivity proposed by Zhang. Experimentally, in 
high-T c cuprates the Zeeman magnetic field can be real- 
ized by applying a magnetic field parallel to ab planes. 
Using non-linear quantum a model and the spherical ap- 
proximation we have found explicit expressions for the 
critical lines and various thermodynamic functions (free 
energy, entropy and specific heat). Obtained results are 
in qualitative agreement with experimentally observed 
features of high-T c cuprates. The maximum of the spe- 
cific heat is depressed and shifted towards lower temper- 



ature in higher magnetic fields. However, it should be 
stressed that the case of the Zeeman field considered in 
the present paper neglects the orbital effects, which are 
important in high-T c cuprates, when the magnetic field 
is applied perpendicularly to ab planes. 
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FIG. 1: T — x — B phase diagram for uj — 3 and w/J — 1. 
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FIG. 2: Temperature dependence of specific heat for various 
magnetic fields, for uj — 3 and w/J = 1. 



